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Abstract 

We model the influence of scattering along the channel and extension regions of dual gate nanotran- 
sistor. It is found that the reduction in drain current due to scattering in the right half of the channel 
is comparable to the reduction in drain current due to scattering in the left half of the channel, when 
the channel length is comparable to the scattering length. This is in contrast to a popular belief that 
scattering in the source end of a nanotransistor is significantly more detrimental to the drive current than 
scattering elsewhere. As the channel length becomes much larger than the scattering length, scattering in 
the drain-end is less detrimental to the drive current than scattering near the source-end of the channel. 
Finally, we show that for nanotransistors, the classical picture of modeling the extension regions as simple 
series resistances is not valid. 
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I. Introduction 

Experimental and theoretical work on nanotransistors has been a hot area of research 
because of significant advance in lithography. The significant advances in lithography have 
led to the construction of nanotransistors with channel lengths smaller than 25 nanometers 
(nm) [1], [2], [3]. It is believed that devices with channel lengths equal to 10 nm may 
become possible in research laboratories [4]. In these nanotransistors, the length scales of 
the channel, gate, screening and scattering lengths, begin to become comparable to one 
another. This is not the case for long channel MOSFETs, where the channel and gate 
lengths are much larger than the scattering lengths. As a result of the comparable length 
scales, it is expected that the physics of nanotransistors will begin deviating from that of 
long channel transistors. 

The resistance of a MOSFET (Fig. 1) with a long channel length can be qualitatively 
thought of as arising in four regions, Extension regions near the source (Ex-s) and drain 
(Ex-d), Channel (Ch), and Contacts. It is believed that the resistance of the contacts 
and extension regions are extrinsic series resistances [5], while the channel resistance is 
intrinsic to the MOSFET. For a given doping distribution, both electrostatics and scat- 
tering of the current carriers play an important role in determining the drive current. 
Electrostatics dictates that the total carrier density in a long channel MOSFET is ap- 
proximately C ox (Vg — Vs), where Vq and Vs are the gate and source voltages. The role 
of scattering of current carriers in long channel transistors is modeled using the mobil- 
ity. For nanotransistors with ultra short channel lengths, there are some deviations in 
the electrostatics from the long channel case [5]. The role of scattering is however not 
well understood in nanotransistors. Most work on nanotransistors use the drift diffusion 
equations which are applicable to long channel MOSFETs or fully ballistic calculations 
based on the Schroedinger equation. A detailed understanding of the influence of scatter- 
ing is important as it is crucial in determining the on-current of nanotransistors. The role 
of scattering is however not straight forward to determine without a calculation because 
scattering tends to change the carrier and current densities in the channel and exten- 
sion regions, both spatially and energetically. Further, the physics of this redistribution 
depends sensitively on the channel and scattering lengths as demonstrated in this paper. 
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The aim of this paper is to model the exact influence of scattering at different spatial 
locations along the channel and extension regions of silicon n-MOSFETs. We consider 
only electron-phonon scattering, which is an important scattering mechanism in devices 
with undoped channels. References [6] and [7] have recently pointed out that electron- 
electron and plasmon scattering may play an important role in degrading nanotransistor 
characteristics. Electron-electron scattering in the drain side will lead to carriers having an 
energy larger than the source injection barrier. The resulting small tail of hot carriers [8] 
will be reflected back into the source-end, there by causing an increase in the source 
injection barrier and a corresponding decrease in drain current. The modeling of these 
effects and interface roughness is beyond the scope of our current work. 

In our calculations, we consider the dual gate MOSFET [9], [10], which is considered to 
be a promising candidate for nanotransistors. The reason for this is the large on-current 
and better scaling properties it offers, when compared to bulk-type MOSFETs [11], [12], 
[13], [14], [15], [16]. 

The outline of the paper is as follows. In section II, we present our simulation results 
on the role of scattering in nanotransistors, where we show that scattering is important 
throughout the device, and not just in the source-end. This is followed by a discussion ex- 
plaining why drain-end scattering is important in nanotransistors (section III). In section 
IV, we show that scattering in the extension regions cannot be modeled as simple series 
resistances. We conclude in section V. All details of our method and approximations are 
given in the appendix. 

II. Where is scattering important?: Simulation Results 

Two devices (Fig. 1) are simulated with the following parameters: 
Device A (similar to the Purdue dual gate MOSFET [17].): Channel length (L Ch ) = 10 nm, 
channel extends from -5 nm to 5 nm, channel thickness (Tch) = 1-5 nm, oxide thickness = 
1.5 nm, gate work function = 4.25 eV, doping in the extension regions = 1 E+20 cm -3 , no 
doping in the channel, drain voltage (Vn) = gate voltage (Vq) = 0.6 V, and the dielectric 
constant of the oxide (e ox )=3.9. 

Device B : Same as Device A, except that L C h = 25 nm, channel extends from - 12.5 nm 
to 12.5 nm and Vg = 0.56 V. In all simulations involving this device, scattering is included 
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only in the channel. The gate length is equal to the channel length for both devices A and 
B. The temperature assumed in all calculations in this paper is 300K. 

We first discuss device A. To elucidate the role of scattering in different spatial regions, 
we calculate the drain current (Id) as a function of the right boundary of scattering, 
YR~scatt- Scattering is included from the edge of the source extension region (-20 nm) to 
YR-scatt in Fig. 2. The ballistic current is 1.92 mA/jum, the value at Yf>_scatt = — 20 nm. 
The channel extends from -5 nm to +5 nm. The main points of this figure are: 

(i) The decrease in current from the ballistic value due to scattering in the source 
extension, channel and drain extension regions are 11.5%, 15.5% and 4% respectively. 
These values point to the well appreciated result that either reducing the length or flaring 
the source extension region will make a nanotransistor significantly more ballistic. 

(ii) The decrease in drain current due to scattering over the entire channel is important. 
That is, scattering in the right half of the channel (0 nm to 5nm) is almost as important 
as scattering in the left half of the channel (-5 nm to nm). 

(iii) The drain current continues to decrease significantly due to scattering in the drain 
extension region. An important question is if this decrease is simply a series resistance 
effect (see section IV). 

We now present results for device B, whose channel length is two and a half times larger 
than device A. The scattering times are nearly the same for the two devices. As a result 
of the larger channel length, the probability for a carrier to energetically relax is larger. 
Here, we find that scattering in the left (-12.5 nm to nm) and right (0 nm to 12.5 nm) 
halves of the channel reduces the drain current by 32% and 15% respectively from the 
ballistic value, and the over all ballisticity (ratio of Current with scattering to Ballistic 
current) is 53% (dashed line of Fig. 4). Again, this points to the importance of scattering 
in the drain-end. 

In lieu of simulating devices with longer channel lengths, we increase the scattering rate 
of device B. The scattering rate is increased by a factor of five by artificially increasing 
the values of the deformation potential quoted in reference [18] by a factor of y/E. Note 
that device B has almost no DIBL and that we self-consistently solve the Green's function 
and Poisson's equations with the larger deformation potentials. The ballisticity of device 
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B with the larger scattering rate is 38%, and the current decreases by 60% and 12% of 
the ballistic value due to scattering in the left and right halves of the channel respectively 
(solid line of Fig. 4). It is also apparent from Fig. 4 that the effect of scattering on drain 
current becomes relatively smaller as Y R scatt approaches the drain-end (12.5 nm). 

III. Discussion 

The results of section II show that scattering at all locations in the channel is important 
in determining the drain current of nanoscale MOSFETs. We first discuss device A. For 
device A, the scattering time (h/2\Im {p r p honon) I) at an energy of E b + 26 meV is 50 fs 
and 24 fs in the source and drain-ends respectively. The scattering times are comparable 
to the semiclassical transit time of 26 fs (Table I). The scattering (11 nm) and channel 
lengths (10 nm) are hence comparable (Table I). It is interesting to note that for this 
device, the argument that the energetic redistribution of electrons in the channel to states 
with kinetic energy in the transport direction well below Ef, will make drain-end scattering 
ineffective fails. 

To understand why drain-end scattering is important for the parameters in device A, 
it is useful to plot the change in barrier height (E b ) with Y R _ Sc att- Fig. 2 shows E b as a 
function of Y R -scatt- It is noted that E b first decreases and then increases, with increase 
in Y R _scatt- The decrease of E b for —20 nm < Y R _ S catt < —4 nm is due to the potential 
drop in the source extension region arising from the increasing series resistance. Note that 
the location of the source injection barrier (Y b ) is -4 nm (Fig. 3). For Y R _scatt > Y b , 
E b increases with Y R scatt- The reason for the increase in E b are the electrons reflected 
towards the source from the right of Y b . Electrostatics, more or less demands that the 
charge in the gate should be approximately C ox (Vg — Vs) [19], [20], like in long channel 
MOSFETs [5]. So, E b floats to higher energies to compensate for the increase in electron 
density from the reflected electrons. This increase in E b contributes significantly to the 
decrease in the drain current even due to scattering in the right half of channel (0 nm to 5 
nm). The increase in E b with increase in Y R scatt becomes smaller in the right end of Fig. 
2 because the electrons scattered here contribute less significantly to the channel charge, 
as will be apparent from the discussion below. 

We now discuss device B. Device B is different from device A in that its channel length is 
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two and a half times longer than that of device A. The importance of scattering in the right 
half of the channel is obvious for device B from the dashed line of Fig. 4. Here, scattering 
in the left (-12.5 nm to nm) and right (0 nm to 12.5 nm) halves of the channel reduce 
the drain current by 32% and 15% respectively from the ballistic value. To complement 
the discussion of device A in terms of we will discuss device B in terms of another 
useful quantity: J(Y, E), which is the current distribution as a function of total energy E 
at Y . J(Y, E) gives us partial information about the energetic redistribution of current 
due to scattering (see end of current section). When the channel length is comparable to 
the scattering length, J(Y, E) is peaked in energy above E b: in the right half of the channel 
(Fig. 5 (a)). Scattering causes reflection of this current towards the source. This is the first 
reason for the reduction in drain current. The second reason is that the reflected electrons 
lead to an increase in the channel electron density (classical MOSFET electrostatics). As 
the charge in the channel should be approximately C ox (Vq — Vs), the source injection 
barrier E b floats to higher energies to compensate for the reflected electrons. The increase 
in Ef, leads to a further decrease in drain current due to scattering in the right half of the 
channel. 

To gain further insight into the role of carrier relaxation, we now discuss device B when 
the scattering length is five times smaller. The scattering length L scaU is defined in Table 
I. Scattering in the right half of the channel for L scaU = 2.2 nm is significantly less 
detrimental to the drain current relative to scattering in the left half of the channel, when 
compared to the device with L scatt = 11 nm. As Lch (25 nm) is much larger than L scatt 
(2.2 nm), multiple scattering events now lead to an energy distribution of current that is 
peaked well below the source injection barrier in the right half of the channel as shown in 
Fig. 5 (b). The first moment of energy (mean) with respect to the current distribution 
function, which is defined by ^^f^'f} , is also shown in Fig. 5. This mean also shows 

' J J dEJ(Y,E) ' ° 

that the carriers relax in a manner akin to bulk MOSFETs as a function of Y in Fig. 5 
(b). Carriers reflected in the right half of the channel can no longer reach Y b due to the 
large barrier to the left, and so contribute less significantly to the charge density. Thus, 
explaining the diminished influence of scattering in the right half of the channel relative 
to the left half of the channel, for devices with the channel length much larger than the 
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scattering length. 

The above discussion would be incomplete without discussing the electrostatic potential 
profiles, with and without scattering. The solid line in Fig. 6 is the electrostatic potential 
in the ballistic limit. Increasing Y^scatt from - 2.5 nm to 2.5 nm causes to increase 
because of carriers reflected towards the source. Further increase in Y^scatt to 7.5 nm 
causes very little increase in Eb because scattering in the right half of the channel is less 
effective in changing the channel electron density. The electrostatic potential changes 
appreciably to the right of Y b due to scattering. It is also interesting to note that the 
electrostatic potential drop for Ynscatt = 7.5 nm is linear to the right of Y& compared to 
the ballistic case because of scattering in the channel. 

We now comment briefly on two issues: 

- The quantity Eb — 2kT that has been discussed before in references [17] and [21]. 

- The influence of elastic scattering without any inelastic scattering. 

For devices A and B, the potential profile in the right half of the channel is well below 
E b — 2kT. Yet, scattering in the right half of the channel is detrimental to drain current, 
relative to scattering in the left half of the channel. The reason for this are the hot electrons 
in the right half of the channel that are reflected to the source-end / Y b . However, if the 
scattering rate in the left half of the channel is large enough to energetically relax the 
electrons to energies comparable to Eb — 2kT, then the scattering of these electrons in 
the drain-end are relatively less detrimental to the reduction in drain current because the 
carriers cannot easily gain an energy of few times the thermal energy. This phenomenon of 
the diminished role of scattering in the channel at the drain-end relative to the source-end 
because of thermalized carriers is seen in Fig. 4 (solid line). 

In the presence of elastic scattering processes such as interface roughness scattering, the 
electron does not loose total energy. However, the kinetic energy in the transport direction 
can diminish at the expense of a corresponding gain in The additional density of 

states for scattering that is available in the drain-end in comparison to the source-end will 
also make drain-end scattering less effective than source-end scattering. While we included 
such process in our calculations, the quantity J(Y, E) captures only the effect of change 
in total energy. A physically motivated study quantifying the relative roles of elastic and 
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IV. Failure of the classic series resistance picture for nanotransistors 

We ask the question if scattering in the extension regions is a simple series resistance. 
The classic series resistance picture [5] relates the current in a device with long extension 
regions to the current in the same device without (or with much smaller) extension regions. 
The relationship is particularly simple for the case where the series resistance in the source 
extension region is negligible [5], 

Ih catt (V D ) ~ I™ scatt {V D - 8V D ) , (1) 

where, /^""(Vd) and l^ scatt {y D — 5Vd)Rd) are the drain currents with and without 
scattering in the drain extension region, at drain biases of Vo and V D — 5Ve> respectively. 
8V D = Ip att (V D )R D , is the electrostatic potential drop in the drain extension region, which 
has a series resistance of Rp- To answer the question on the appropriateness of the classic 
series resistance picture, we consider a case where the channel and source extension region 
are ballistic. Scattering is introduced only in the drain extension region with deformation 
potentials that are ^/E times larger than in silicon (scattering time is 5 times smaller). 

Fig. 7 shows the decrease in drain current with Y R ^ S catt- The striking point of Fig. 7 
is the super-linear decrease of drain current. The Id{Vd) curves (inset of Fig. 7) predict 
a significantly smaller decrease in drain current with increase in Y^scatt when Eq. 1 
is used. It is helpful to estimate the drain current from Eq. 1 and compare it to the 
calculated value. For Device A in Fig. 7, the voltage drop in the drain extension region 
with scattering is approximately 100 mV (plot not shown). Now, if Eq. 1 is used to 
estimate the drain current with scattering in the drain extension region and if we take 
SVd = 200 mV, which is larger than the estimated 100 mV, then we find the drain current 
to be 1.83 mA//im (inset of Fig. 7). The calculated drain current is however much lower 
at 1.38 inA/jum! 

The physics of the large reduction in drain current for the smaller values of Y^scatt is 
essentially that discussed in section III: When scattering in the channel does not effectively 
thermalize carriers, the current distribution is peaked at energies above Eb, upon carriers 
exiting the channel. Scattering in the drain extension region then causes reflection of 
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electrons towards the source-end. As a result, increases so as to keep the electron density 
in the channel approximately C 0X (Vq — Vg). The drain current decreases dramatically as 
a result of the increase in E\>. Admittedly, this argument in terms of C ox (Vg — Vs) is over 
simplified but it seems to capture the essential point. The main point is that if carriers are 
not relaxed upon exiting the channel (as would be the case for nano-transistors) , then, the 
drain extension region cannot be modeled by a simple series resistance. That is, Eq. (1) 
fails for nano-transistors where the channel length is comparable to the scattering length. 
The effect of the drain extension region in causing a reduction in drain current would be 
small in the following cases: 

(i) The channel is much longer than the scattering length such that the carriers exiting 
the channel at the drain-end are energetically relaxed / thermalized. Then, the modeling 
of the drain extension region as a simple series resistance would be appropriate. This is 
seen in the right end of Fig. 7, where, upon sufficient relaxation of electrons, the decrease 
in current with increase in Y R ^ Scatt becomes much smaller. 

(ii) The drain extension region rapidly flares out. Then, the probability for a scattered 
electron to return to the source-end will be small due to the larger number of modes 
available in the drain extension region. A careful analysis on how fast the drain extension 
region flares out should also take into account the role of the Miller effect. 

V. Conclusions 

In conclusion, we find that the potential profile, channel and scattering length scales 
play an important role in determining the relative importance of scattering at different 
locations along the channel of a nanotransistor. In devices where the channel length is 
comparable to the scattering length, the role of scattering in the drain-end (right half of the 
channel) is comparable to the role of scattering in the source-end (left half of the channel), 
in reducing the drain current (Fig. 2 and dashed line of Fig. 4). This is contrary to a belief 
that scattering is significantly more important in the source-end of the device. The reason 
for the detrimental role of scattering in the drain-end are the hot carriers in the drain-end. 
When the channel length is much larger than the scattering length, then scattering in the 
source-end becomes relatively more important than scattering in the drain-end (solid line 
of Fig. 4). In this case, we stress that it is the energetic redistribution of carriers due to 
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scattering in the source-end to energies below the source injection barrier (E b ) that makes 
scattering in the drain-end relatively less detrimental to the drain current. 

The classical series resistance picture for modeling the narrow extension regions fail 
for nanotransistors. The reason for this failure are the hot carriers entering the drain 
extension region. A straight forward option to enable the usage of the series resistance 
picture is to push the region treated as a drain series resistance further to the right, such 
that all carriers entering this region are energetically relaxed. A more interesting option 
of altering the classical series resistance picture to account for the hot carriers in the drain 
end of nanotransistors was not considered in this paper. 

The relative importance of scattering in the drain-end of nanotransisors, where the chan- 
nel length is comparable or smaller than the scattering length, points to the importance of 
making the extension regions small. Long extension regions in nanotransistors will affect 
the performance (drive current) much more adversely than in long channel transistors. 

Appendix 

The approach consists of solving the nonequilibrium Green's function and Poisson's 
equations. The effective mass Hamiltonian considered is, 



, 2 



2 r 



d ( 1 d \ d ( 1 d \ d ( 1 cT 

dx \m b x dx J dy \m b dy J dz \m b dz 



+ V(x,y), (2) 



y 

where (m^, m b y , m b ) are the (x, y, z) components of the effective mass in valley b of silicon, 
and the potential does not vary in the z direction. The gate oxides are treated as hard 
walls, the channel is extremely narrow (1.5 nm), the drain and gate biases are smaller than 
0.7 V, and the dual gate FET is perfectly symmetric in the X-direction of Fig. I. The 
first three subband energy levels in the source extension region are approximately equal 
to 173 meV, 691 meV (both due to m y = 0.98m ) and 891 meV (due to m y = 0.19m ) 
above the bulk conduction band. The Fermi energy of bulk silicon at the doping density 
considered (1E+20 cm -3 ) is approximately 60 meV above the conduction band. For the 
doping density considered, electrons are primarily injected from the source into the first 
subband. At the drain end, more than one energy level can in principle contribute to 
current. As only a few subbands are populated, we model transport in these subbands in 
an approximate way using the ID Schroedinger equation as outlined below. We find the 
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spatially dependent subband energies E n (y) by solving Schroedinger's equation at each 
y-cross section (y is only a parameter), 

h 2 d ( 1 d 



~K^Z )+V(x,y) 



V n (x,y) = E n (y)V n (x,y) . (3) 



1m\ dx \m b x dx / 

n = u,b, where v and b represent the quantum number due to quantization in the X- 
direction and the valley respectively. The valley indices b are required in the calculations 
of the self-energies for scattering as will be discussed below. In our calculation, we typically 
retain only the three lowest energy levels. Coupling between the subbands is neglected 
except via phonon coupling. For the device dimensions and voltages considered, reference 
[22] found the approximation of considering decoupled subbands to hold good for ultra 
thin body phase coherent MOSFETs. We solve the following equations for the Green's 
functions, 



h 2 k 2 ( h 2 d ( 1 d \ N 

2m« ^ 2 dy \m»dy) + n{V) , 



G r n (y,y',k z ,E) 

J d Vl E r n (y, Vl , k z , E)G r n (y u y', k z , E) = S(y - y') , and (4) 



~ 2^™ ~ \ Ydy [m^dyj + 



G a n (y,y',k z ,E) 

-/ d Vl i: r n (y,y 1 ,k z ,E)G':(y 1 ,y',k z ,E) = J dy E"(y, y\, k z , E)G^(y 1 , y', k z , E) ,(5) 



where, a G>, <. to™ and m™ are the effective masses of silicon in the y and z directions 
that give rise to subband index n. 

The self-energies, 'E% >,< can be written as, 

K = K,C + K,Phonon' where ( 6 ) 

ya _ ya , ya /y\ 

n, Phonon n,el ~ n,inel " V / 

E" c is the self-energy due to the leads. The phonon self-energy E" Phonon consists of two 
terms, E° el due to elastic and E" inei due to inelastic scattering. The self-energy due to 
the leads is non zero only at the first (source) and last (drain) grid points because gate 
tunneling is neglected. 

The following common approximations to calculate the phonon self-energies are used: (i) 
Phonon scattering is treated only within the self-consistent Born approximation, (ii) The 
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phonon bath is assumed to always be in equilibrium, and so their occupation numbers are 
given by the Bose- Einstein distribution function with a spatially independent temperature, 
(iii) The correlation between subbands n and n! n) are neglected, (iv) Scattering due to 
phonons is assumed to be isotropic. That is, the scattering rate from (k Z) E) to (k' z , E') does 
not depend on k z and k' z . This approximation is computationally advantageous because 
the self-energies due to phonon scattering appear only as diagonal terms in Eqs. 4 and 
5. One can derive from these assumptions that the self-energies due to electron-phonon 
scattering at grid point yi are given by [23], [24], 



\ 171 f 1 



Xf nel , n (yi,E) = Y.D^X-JL / dEz — 

[riB^uj^G^iyi, E z , E - hw v ) + (n B (hu> v ) + l)G>,{y u E z , E + hu> v )] , (9) 



and 



[n B (^)G>,(j/i, E z , E + hu v ) + (n B (hu; v ) + l)G>,(y i: E z , E - hoo v )} . (10) 

a G>, <, r in Eq. 8, r] represents the phonon modes, and the square of the matrix elements 
for phonon scattering are given by, 

1 D 2 kT 

Kn> = + (11) 



DLH . (i2) 



2pu m ' puj fv 

The contribution to elastic scattering is only from acoustic phonon scattering. The values 
of the deformation potential, D A , D gv and Df v , and phonon frequencies u gv and Uf v are 
taken from [18]. p is the mass density, k is the Boltzmann constant, T is the temperature 
and v is the velocity of sound, b and b' are indices representing the valley. The following 
scattering processes are included: acoustic phonon scattering in the elastic approximation 
and g-type intervalley scattering with phonon energies of 12, 19 and 62 meV. It was verified 
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that f-type (19, 47 and 59 meV phonon) intervalley scattering did not significantly change 
our results and conclusions. This can be rationalized by noting that f-type intervalley 
scattering involves subbands with energies higher than the lowest subband. In the regions, 
where scattering was not included, the deformation potential was set equal to zero. 
^inei n can be obtained using the Kramers- Kronig relationship, 

Re 



^neUv.E)] = -P / dE> L ™ h f E ^and (13) 



lm 



1 



AneiM E )\ = Yi F'Wf" E) - E< ei)n (y,, E)\ , (14) 

where P stands for the principal part of the integral. Note that the self-energies due to 
electron-phonon scattering depend only on the total energy E (and not on k z ) due to the 
assumption of isotropic scattering. 

The self-energy due to phonon scattering, has real and imaginary parts, both of which 
vary with energy. The imaginary part of the electron-phonon self-energy which is central 
to our calculations is responsible for scattering induced broadening of energy levels and 
energetic redistribution of carriers. The real part of the self-energy which contributes to 
the shift of the quasi-particle energy levels, appears as a real potential (like the electrostatic 
potential) in the Green's function equations (Eqs. 4 and 5). To evaluate the importance of 
the real part of the self-energy in our calculations, we performed simulations with acoustic 
phonon scattering in silicon, with and without the real part of the self-energy included. 
We find that the drive current calculated with the real part of the self-energy set to zero 
in general agrees to within 2 percent of the current calculated with the real part of the 
self-energy included. This result is not totally surprising because MOSFET electrostatics 
tends to shift the potential profile appropriately to determine the correct charge under the 
gate. In the calculations presented in this paper, the real part of the self-energy is set to 
zero. 

In the numerical solution, we consider N uniformly spaced grid points in the redirection 
with the grid spacing equal to Ay. The discretized form of Eqs. 4 and 5 are: 

A iti G r n ( yi , yl k z , E) + A i>i+1 G r n (y i+1 , & k z , E) + A^G^y^, y[, k z , E) = |^ , and (15) 
Ai,iG°(yi, y[, k g , E) + A iti+1 G%(y i+1 , y[, k z , E) + A^G^y^, y[, k z , E) = 
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^( yi ,E)G a n (yi,y'i,k z ,E) , (16) 

where, 

A hl = E-^--^- 2 -E n ( yi )-^ n (y l: k z ,E)^d (17) 

z y <j 

= (18) 

The self-energy due to the source and drain leads contribute only to grid point 1 (left 
end of the source extension region) and grid point N (right end of the drain extension 
region), and are given by [25]: YT nC (yi, k z , E) = ( 2m f Ay2 ) 2 g s (k z ,E), Y7 nC (y N , k z , E) = 
(2^) 2 9d(k z ,E), E< c ( yi ,k z ,E) = -2iIm(££ jC (yi, k z , E))f s (E), E< c (y N ,k z ,E) = 
-2ilm(Z^ c (y N , k z , E))f d (E), Z> c (yi, k z , E) = 2iIm(E^ c ( J/l , k z , E))[l - f s (E)}, and 
^n,c(VN,k z ,E) = 2ilm(T, r nC (y N: k z ,E))[l - f d (E)], where y 1 an y N are the left (source- 
end) and right (drain-end) most grid points respectively, g s (k z ,E) and g d (k z ,E) are the 
surface Green's functions of the source and drain leads respectively, and f s and f d are the 
Fermi functions in the source and drain contacts respectively 

The non equilibrium electron and current densities are calculated in both the channel 
and extension regions using the algorithm for G K in [26], which avoids full inversion of 
the A matrix. For completeness, we state the expressions for the electron and current 
densities used [26], 

n n ( yi ,k z ,E) = -iG<(yi,y h k z ,E) (19) 
e h 2 

J n { yi ,k z ,E) = - 2m n Ay 2 \ G niVh Vi+i, k z , E) - G<(y i+1 , Vi , k z , E)] . (20) 

Note that Eqs. 19 and 20 do not include spin and valley degenaracies. The total electron 
and current densities at grid point yi are given by, 



= %^lfJ iE 'iw:^ E - E) (21) 



rhV2 

™2 rdE r l 



JM = %i^!i;! iE '7K Uy '' E '' E) ' (22) 

where the prefactor of 2 in the above equations account for two fold spin degenaracy. 
While the transport equations are solved in one dimension, we solve Poisson's equation 
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in two dimensions. The two dimensional electron density used in Poisson's equation is 
computed from Eqs. (3) and (19) using, 



The boundary conditions to Poisson's and Green's function equations are applied at the 
ends of the source and drain extension regions (left and right ends of the source and drain 
extension regions shown in Fig. 1). In solving the the Green's function and Poisson's 
equation, note that an applied bias corresponds to a difference in the Fermi levels used in 
the source and drain regions. The electrostatic potential at the left and right most grid 
points of the source and drain extension regions respectively are calculated self consistently 
using the boundary conditions. 

Finally, we make a comment on the need for solving quantum mechanical equations 
to capture the essential effect of hot carriers, described in this paper. The phase of the 
electron is not central to the physics described in our paper (though the exact value of the 
drain current depends on it). In calculating the drain current, the quantum mechanical 
effects of quantization in the X-direction and tunneling along the Y-direction (Fig. 1) 
can be accounted for semiclassically. So, we feel that a method such as the Monte Carlo 
method approach to nanotransistors [27], which keeps track of the details of the energetic 
redistribution of electrons at various spatial locations, will well describe many aspects of 
the role of scattering. 
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Fig. 1: Schematic of a Dual Gate MOSFET (DG MOSFET). Ex-s and Ex-d are the 
extension regions and the hatched region is the channel. The white region between the 
source / drain / channel and gate is the oxide. The device dimension normal to the page 
is infinite in extent. 

Fig. 2: Plot of drain current (Id) versus the right boundary of scattering (Y R _scatt) 
for device A. The scattering time is comparable to the transit time through the channel. 
Scattering is included from -20 nm to Y R _s C att- Note that scattering in the right half of the 
channel (0 nm to 5 nm), which is to the right of the '/c^T layer', is almost as deleterious 
to current flow as scattering in the left half of the channel (-5 nm to Onm). The black 
crosses represent E b as a function of Y R _ Sc att- Inset: Ballistic I D versus V D for V G = 0.6 
V, showing substantial DIBL. Scattering is included both in the channel and extension 
regions. 

Fig. 3: Energy of the lowest subband (E ± ) versus Y for device A in the ballistic limit. 
Eb and Y b are the energy and position of the source injection barrier respectively. Potential 
= -Ml 

e 

Fig. 4: Plot of drain current versus Y R _ Sc att for device B. Scattering is included from 
-12.5 nm to Y R ^scatt- For L scatt = 11 nm (dashed line) and 2.2 nm (solid line), the effect 
of scattering in the right half of the channel (0 nm to 12.5 nm) corresponds to nearly a 
third and sixth respectively of the total reduction in drain current. This figure points to 
the relatively smaller role of drain-end scattering in comparison to source-end scattering, 
when Lch becomes much larger than L scaU - Scattering is included only in the channel for 
both cases. 
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Fig. 5: The solid lines represent J(Y,E) for Y equal to -17.5, -12.5, -7.5, -2.5, 2.5, 
7.5, 12.5 and 17.5 nm, from left to right respectively. The dashed lines represent the first 
resonant level (Ei) along the channel. The dotted lines represent the first moment of energy 
(mean) with respect to the current distribution function J(Y, E), which is j dEJ (YE) • ( a ) 
and (b) correspond to L scatt = 11 and 2.2 nm respectively in device B. Scattering is 
included every where in the channel, (a) and (b) correspond to the Ynscatt — 12.5 nm 
data points of the dashed and solid lines of Fig. 4. Scattering is included everywhere in 
the channel but not in the extension regions. 



Fig. 6: Electrostatic potential versus Y for device B. Scattering from -12.5 nm to 2.5 
nm causes a large change in the source injection barrier (E b ). Scattering to the right of 
2.5nm causes a much smaller change in E b . In the absence of scattering, the potential 
profile in the channel tends to flatten. The potential drop (or E\) along the channel is 
more ohmic / linear in the presence of scattering. 



Fig. 7: Id versus Y R _ Scatt for device A with scattering present only in the drain 
extension region from 5 nm to 30 nm. The large reduction in drain current is due to 
scattering of hot carriers from the drain extension region back in to the channel. The 
physics of this effect is completely different from 'classical series resistance' in MOSFETs, 
which is a much smaller effect. The results obtained from the 'series resistance' and 
'scattering calculations' (this paper) are indicated by the arrows. The electron-phonon 
scattering time is five times larger than in Fig. 2. Inset: Drain current versus drain 
voltage in the ballistic limit, showing the drain current estimate from the series resistance 
picture and from our calculation. 



February 1, 2008 



DRAFT 



IEEE TRANSACTIONS IN ELECTRON DEVICES 



20 





Device A 


Devic B 


T scatt at source-end (s) 


5.0 E-14 


5. 0E- 14 (1.0 E-14) 


x scatt at drain-end (s) 


2.5 E-14 


2.4 E-14 (4.8 E-15) 


T transit at E b +26 meV (s) 


2.6 E-14 


6.4 E-14 


transit at 60 meV (s) 


2.0 E-14 


5.6 E-14 


v at Y b , E b +26 meV (m/s) 


2.2 E+5 


2.2 E+5 


v at Y b , 60 meV (m/s) 


3.5 E+5 


2.8 E+5 


L Ch (nm) 


10 


25 


v*(Xcatt at Y b ) 
E=E b +26meV (nm) 


11 


11 (2.2) 


• x scatt - scattering time (hbar/2Im(E r )) 

• Vansit " shortest semiclassical transit time for electron with a given total energy = integeral 
jdy/[2(E-V(y))/m] 1/2 

• v - semiclassical velocity at y = [2(E-V(y))/m] 

• For Device B, quantities in brackets are for the case of five times larger scattering rate 



TABLE I 

Estimates of scattering time, transit time, velocity and scattering length. 
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Fig. 2 
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Fig. 3 
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Device A, L_. = 10 nm, L = 2.2 nm 




Fig. 7 
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